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2 Abstract Zhou and Cess [2001] developed an algorithm for retrieving surface 

3 downwelling longwave radiation (SDLW) based upon detailed studies using radiative 

4 transfer model calculations and surface radiometric measurements. Their algorithm 

5 linked clear sky SDLW with surface upwelling longwave flux and column precipitable 

6 water vapor. For cloudy sky cases, they used cloud liquid water path as an additional 

7 parameter to account for the effects of clouds. Despite the simplicity of their algorithm, it 

8 performed very well for most geographical regions except for those regions where the 

9 atmospheric conditions near the surface tend to be extremely cold and dry. Systematic 

10 errors were also found for scenes that were covered with ice clouds. An improved version 

11 of the algorithm prevents the large errors in the SDLW at low water vapor amounts by 

12 taking into account that under such conditions the SDLW and water vapor amount are 

13 nearly linear in their relationship. The new algorithm also utilizes cloud fraction and 

14 cloud liquid and ice water paths available from the Cloud and the Earth’s Radiant Energy 

15 System (CERES) single scanner footprint (SSF) product to separately compute the clear 

16 and cloudy portions of the fluxes. The new algorithm has been validated against surface 

17 measurements at 29 stations around the globe for Terra and Aqua satellites. The results 

18 show significant improvement over the original version. The revised Zhou-Cess 

19 algorithm is also slightly better or comparable to more sophisticated algorithms currently 

20 implemented in the CERES processing and will be incorporated as one of the CERES 

21 empirical surface radiation algorithms. 

22 
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1 1. Introduction 


2 Using detailed studies based upon radiative transfer model calculations and surface 

3 radiometric measurements, Zhou and Cess [2001] formulated algorithm development 

4 strategies for retrieving surface downwelling longwave radiation (SDLW). Their studies 

5 demonstrated that clear sky SDLW could be largely determined by surface upwelling 

6 longwave flux and column precipitable water vapor. For cloudy sky cases, they used 

7 cloud liquid water path as an additional parameter to account for the effects of clouds 

8 instead of using cloud base height which is theoretically a more direct factor in 

9 determining the SDLW. An illustrative algorithm was derived and tested using 

10 observational data from Atmospheric Radiation Measurements (ARM) Program [Stokes 

11 and Schwartz 1994] measurements at the U.S. Southern Great Plain (SGP) and Tropical 

12 Western Pacific (TWP) sites. 

13 Since the algorithm was derived and tested for mid-latitude and tropical conditions, 

14 there was concern that the algorithm might not perform well for extremely cold and dry 

15 conditions. Indeed, once data had become available from the ARM North Slope of 

16 Alaska (NSA) [Stamnes et al., 1999], large biases were found when the algorithm was 

17 applied to surface measurements from that location. Sensitivity studies demonstrated that 

18 the algorithm significantly underestimates SDLW when atmospheric water vapor was 

19 low. 

20 Meanwhile, the Clouds and the Earth’s Radiant Energy System [CERES; Wielicki et 

21 al. 1996] Surface Radiation Budget (SRB) team started to test the algorithm for possible 

22 global application. The CERES program is designed to provide crucial cloud and 

23 radiation measurements for studying cloud-radiation interaction. CERES instruments 
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were launched aboard the Tropical Rainfall Measuring Mission (TRMM) in November 

2 1997 and on the Earth Observation System (EOS) Terra satellite in December 1999 and 

3 Aqua satellite in 2002. The space-borne CERES radiometers provide broadband total, 

4 shortwave (SW) and infrared window measurements at the top of atmosphere (TOA). 

5 Differencing the SW from the total measurements allows for a derivation of the longwave 

6 (LW) value. Deriving reliable estimates of SRB parameters is another important objective 

7 of the CERES project whose goal is to provide a complete picture of the energy budget of 

8 the earth-atmosphere system. Since the SRB cannot be directly measured by satellite- 

9 borne instruments, the surface fluxes are derived with several different methods using 

10 combinations of radiation models, data assimilation products, and satellite measurements. 

11 The Surface and Atmospheric Radiation Budget [SARB; Charlock et al. 1997] 

12 component of CERES represents one such method where shortwave and longwave fluxes 

13 at the surface, at three levels in the atmosphere, and at the TOA are computed with a 

14 radiative transfer model. In addition to SARB, surface fluxes are being derived within 

15 CERES using two SW and two LW models, which are based on TOA-to-surface transfer 

16 algorithms or fast radiation parameterizations. These models are the Li et al. [1993] 

17 model (SW model A, clear sky only), the Darnell et al. [1992] model (SW model B), the 

18 Inamdar and Ramanathan [1997] model (LW model A, clear sky only), and the Gupta et 

19 al. [1992] model (LW model B). These models were incorporated into CERES products 

20 to provide independent sources of surface fluxes to compare with SARB results [Gupta et 

21 al. 2004]. The Zhou-Cess algorithm represents a new methodology for deriving SDLW 

22 globally for both clear sky and cloudy sky using parameters readily available from 

23 satellite measurements. A vigorous test was performed on instantaneous flux derived by 
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1 Zhou-Cess algorithm as part of single-scanner footprints (SSF) products on board Terra 

2 and Aqua satellite against matched ground measurements from surface radiation 

3 measurement sites around the globe. Despite the simplicity of their algorithm, it 

4 performed very well for most of the geographical regions. Large biases, however, were 

5 found for certain regions, most notably the Polar Regions, where the atmosphere is 

6 extremely cold and dry. In addition, systematic errors were found for regions covered 

7 with ice clouds. 

8 The modifications to the algorithm discussed in the present work are aimed at 

9 addressing the low water vapor and ice cloud situations with the goal of making the 

10 algorithm applicable for global implementation in the CERES processing. In section 2, a 

1 1 brief description of the original algorithm will be given followed by bias analysis of the 

12 algorithm when applied to ARM NS A data. In section 3, the algorithm is revised using 

13 collocated CERES cloud parameters and surface radiation measurements. The new 

14 algorithm is validated for CERES Terra and Aqua satellites with surface measurements 

15 from global network. A summary of this work is given in section 4. 

16 

17 2. Low Water Vapor Bias 

18 2.1 Original algorithm 

19 The Zhou-Cess algorithm took the form: 

20 SDLW = a + b • SULW + c • ln(PWV) + d • [In(PWV)] 2 + e • ln(l + f • LWP) (1) 

21 where SULW is the surface upwelling longwave flux computed from the 2-meter air 

22 temperature using Stefan-Boltzmann’s law assuming an emissivity equal to unity. PWV 

23 is the column precipitable water vapor and LWP is the cloud liquid water path, both in 


5 



1 centimeters. The regression coefficients a, b c, d, e, f have values of 123.86, 0.444, 56.16, 

2 -3.65, 5.30, 1226.0, respectively. The algorithm was formulated based on detailed studies 

3 of radiative transfer models and observational data. The unique feature of this algorithm 

4 is that it uses one formula to compute the clear sky and all sky SDLW. The LWP is used 

5 as a surrogate of cloud base height since LWP is found to be correlated with cloud base 

6 height based upon observational data and upon thermodynamic arguments by Hack 

7 [1998]. LWP is also better defined for any spatial and temporal grid. The algorithm was 

8 derived using observational data obtained from six Intensive Observation Periods (IOP) 

9 at the ARM SGP site and was verified with nine other SGP IOP datasets and TWP Manus 

10 data [Zhou and Cess 2001], The algorithm was not tested with observational data from 

1 1 other geophysical regions because of data availability problems. 

12 

13 2.2 The ARM NS A data 

14 Recent climate modeling and diagnostic studies indicate the Polar Regions are 

15 particularly sensitive to global climate change and important to mid-latitude climate and 

16 weather systems [Stamnes et ah, 1999]. Radiation tends to dominate the Arctic heat 

17 budget in all seasons. Due to the extreme weather conditions experienced at high 

1 8 latitudes, the algorithms developed for mid-latitudes frequently do not work well for the 

19 Polar Regions. Thus, prudence necessitated testing the Zhou-Cess algorithm once high 

20 quality observations from ARM NSA became available. The data used in the study was 

2 1 from January to December, 2000 from the ARM Barrow facility, which is located at the 

22 northernmost point (71.32°N, 156.61°W) in the United States, 528 km north of the Arctic 

23 Circle. The data were taken from the same instruments and processed in the same manner 
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1 as those used by Zhou and Cess [2001] for the SGP and TWP sites. The SDLW fluxes 

2 were pygeometer measurements and the surface upwelling fluxes were computed from 

3 the Surface Meteorological Observation Station (SMOS) 2-m surface air temperature. 

4 The column precipitable water and cloud liquid water were both measured by microwave 

5 radiometers (MWR). All data were averaged into half-hour products. The results show 

6 that the algorithm mostly underestimates the SDLW, with large negative bias at the low 

7 SDLW (Fig. 1). Further analysis found that 94% of the underestimated cases are related 

8 to very low water vapor amount (PWV < 0.81 cm). The low water vapor amount has 

9 resulted in very large negative value in the ln(PWV) term since the logarithmic function 

10 decreases very rapidly with decreasing water vapor below 1 cm. To remedy the problem, 

1 1 the ln(PWV) tenn was replaced by ln(l+PWV) to approximate a near-linear relationship 

12 between PWV and SDLW when water vapor amounts are very low. This modification 

13 has the added advantage of producing a very simple and reasonable relationship between 

14 the SDLW and the water vapor term which prevents the PWV from producing a negative 

15 contribution to the flux as the water vapor amounts asymptote toward zero. 

16 Although Zhou and Cess [2001] conducted radiative transfer calculations for 6 

17 default Moderate -Resolution Transmittance Radiation Model (MODTRAN, Wang et al. 

18 1996) atmospheres (including subarctic summer and subarctic winter atmosphere) and 

19 their variations for providing basic relationships between SDLW and other parameters, 

20 the actual algorithm was derived using only observations at SGP site for illustrative 

21 purpose. For any nonlinear relationship, accounting for the full range of conditions is 

22 critical to derive statistical relationships that apply to most situations. This is because the 

23 statistical relations (usually derived with least square fitting) will lean toward highly 
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1 sampled situations and miss the under-sampled situations. 

2 The importance of sampling can be illustrated in Fig. 2 where a linear relationship is 

3 calculated for MODTRAN computed SDLW and SULW from Zhou and Cess [2001] 

4 (upper panel). Using the upper portion (PWV>=0.5 cm) and the lower portion (PWV < 

5 0.5 cm) of the data gives very different slope and offset for a linear relationship. The log- 

6 square fit of PWV to the ratio of fluxes also depends significantly on different portions of 

7 the sampling data (lower panel). The reason why the SGP algorithm does not produce 

8 good results for NSA data is because SGP data only represent middle range of the data, 

9 and thus, is not applicable to the low end of the curve. Fig. 3 illustrates the probability 

10 distribution functions (PDF) of SULW, SDLW and PWV from global networks of 29 

1 1 sites (Fig. 4) and those from ARM SGP and NSA sites. An examination of Fig. 3 reveals 

12 that the dynamical range of surface longwave fluxes and PWV from SGP miss a small 

13 fraction of data where SULW is larger than 500 Wm' 2 or PWV greater than 5 cm as 

14 shown in the global PDFs. More importantly, however, the NSA PDFs clearly 

15 demonstrate that the distributions of these quantities largely misrepresent regions of 

16 extremely cold and dry areas. The misrepresentation of low water vapor situations is 

17 particularly serious because of nonlinear logarithmic relationship between water vapor 

18 and SDLW. 

19 Due to data availability, the original algorithm only used cloud liquid water path to 

20 account for cloud effects to the SDLW. For the SGP and TWP site, the effect of ice 

21 clouds on the SDLW is relatively small since most of the ice clouds are at high altitudes. 

22 Such is not the case for the polar region, where most of the clouds are ice cloud located at 
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1 low altitudes. Moreover, since the atmospheric water vapor concentrations for the Polar 

2 Regions are also very low, the effect of ice clouds on SDLW is not negligible. 

3 

4 3. Satellite Implementation 

5 3.1 Data 

6 The CERES Single Scanner Footprint (SSF) product contains one hour of 

7 instantaneous CERES data for a single scanner instrument. The SSF product combines 

8 instantaneous CERES data with scene information and cloud properties derived from a 

9 high spatial resolution imager such as the Visible/Infrared Scanner (VIRS) on TRMM or 

10 the Moderate-Resolution Imaging Spectroradiometer (MODIS) on Terra and Aqua 

11 satellites. The cloud properties are computed in the cloud subsystem of CERES 

12 processing [Minnis et al. 1997]. Thus, all of the input parameters for the Zhou-Cess 

13 algorithm are already computed or assembled in the current SSF processing. 

14 The ground measurements are taken from the CERES/ARM Validation Experiment 

15 [CAVE; Rutan et al. 2001] database which is maintained at NASA LaRC in a Web- 

16 accessible fonn for use in the CERES project and is available to the outside science 

17 community. The surface observations used in this study include 5 ARM/SGP sites, 1 

18 ARM/TWP site, 7 National Surface Radiation Budget Network (SURFRAD; Augustine 

19 et al., 2000) sites, 8 Baseline Surface Radiation Network [BSRN; Ohmura et al. 1998] 

20 sites, 6 Climate Monitoring and Diagnostics Laboratory (CMDL) sites, the LARC/COVE 

21 site [Jin et al., 2002], and one NREL [Myers et al. 1999] site (Fig. 4). These sites were 

22 selected based on data availability and statistical representative of different geographical 

23 regions, i.e., ocean, land, desert, arctic regions. Temporal matching of the satellite and 
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1 site fluxes was done at the highest resolution of the site data. Spatial matching was done 

2 to a distance of 10 km between the location of the site and the center of the CERES 

3 footprint. Values for all CERES footprints within the 10 km range of the sites and within 

4 the 1 minute interval were averaged together for comparison with the corresponding 

5 ground-based values. 

6 For the purpose of deriving, as well as validating, a new algorithm, it is necessary to 

7 separate the training and validating data set. Here we have chosen 1 5 ground stations and 

8 collocated cloud parameters from Terra satellite during 58 months from March 2000 to 

9 Dec 2004 as a training data set. All 29 ground stations and their collocated Terra products 

10 for the same period are used in the validation since only clear and overcast portions of the 

1 1 training data set are used in the derivation of the algorithm. The collocated Aqua and 

12 ground measurements from all 29 stations covering the period from July 2002 to March 

13 2005 serve as additional data set for validation. 

14 

15 3.2 Performance of the original algorithm 

16 Fig. 5 shows the scatter plots of SDLW computed by the original Zhou-Cess 

17 algorithm for Terra satellite and collocated ground measurements from all 29 sites. For 

18 most geographical regions (continental, desert, island and coastal area), the algorithm 

19 performs reasonably well; however, there is a significant underestimation of the SDLW 

20 in the Ant( arctic), which is defined to include both the Arctic and Antarctic, when SDLW 

21 is below 200 Win'”. Most of the low values were observed at the South Pole where mean 

22 water vapor is only 0.34 cm. When the data are stratified with clear (f c i r > 0.999), water 

23 cloud (LWP > 5 g/m 2 ) and ice cloud (LWP < 5 g/m 2 and ice water path (IWP) > 1 g/m 2 ) 
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1 cases, we find that there is positive bias for most of the clear sky cases except in some 

2 Ant(arctic) cases (Fig. 6). Large negative bias was found for both water cloud and ice 

3 cloud, with ice clouds having the largest negative bias. The large negative bias for cloudy 

4 sky flux might be due to lower cloud liquid water path generally observed from satellite 

5 than those observed from MWR at the SGP site. The PDFs of Fig. 7 show that the 

6 satellite rarely observes LWP greater than 1 mm over the SGP site while MWR observe 

7 many cases with LWP greater than 1 mm. We also found that some of the MWR 

8 measured cloud liquid water used in deriving the original algorithm might have been 

9 contaminated by rain or wetness from other forms of precipitation due to their unlikely 

10 larger value (LWP >1.3 mm). The systematic difference between ice cloud and water 

1 1 cloud also indicates that ice clouds can play an important role for conditions other than 

12 the wann, moist tropical conditions. Taken together, these observations suggest that the 

13 algorithm should be re-derived using satellite observed cloud parameters and using more 

14 sites globally so that the algorithm better represents of various atmospheric conditions. In 

1 5 the revised algorithm, the effect of ice cloud will be considered by including ice water 

16 content in the same manner as cloud liquid water although with a smaller weight. 

17 

18 3.3. Revised algorithm 

19 Based on the above analysis, a new algorithm was formulated which computes the 

20 fluxes for clear portions (F c i r ) and cloudy portions (F c id) of the sky separately and then 

2 1 sums the results for the all sky flux (F a u): 

22 F dr = aO + al • SULW + a2 • ln(l + PWV) + a3 • [ln(l + PWV)] 2 (2) 

23 F dd = bO + bl • SULW + b2 • ln(l + PWV) + b3 • [ln(l + PWV)] 2 + 
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b4 • ln(l +LWP) + b5 • ln(l+IWP) 


( 3 ) 


2 Fa// - F c l r • fd r + F tfd • (1.0 - f c lr) (4) 

3 F net = F aU - SUL W (5) 

4 aO = 37.687, al = 0.474, a2 = 94.190, a3 = -4.935 

5 bO = 60.349, bl =0.480, b2 = 127.956, 

6 b3 = -29.794, b4= 1 .626, b5 = 0.535 

7 where f c i r is the fraction of clear area in a single CERES footprint. SULW and PWV 

8 follow the same units as in (1). LWP and IWP (in g/m") are not representative of the total 

9 for the pixel, but only for the cloudy portion, since the cloud fraction is explicitly taken 

10 into consideration. The cloud fraction and cloud water path are part of SSF product and 

1 1 derived from CERES cloud analysis. The column precipitable water vapor is integrated 

12 from CERES meteorological data. Moreover, clear sky is defined when f c i r is greater than 

13 0.999 and both LWP and IWP are set to zero. SULW is computed using a surface skin 

14 temperature with an emissivity of unity. Test cases run to incorporate a global emissivity 

1 5 map demonstrated no improvement in the algorithm. 

16 As noted previously, the above regressions are derived using satellite observations 

17 from the Terra platform and matched ground measurements from 15 sites around the 

18 globe (Fig. 4) from March 2000 to Dec 2004. There are 6028 clear sky cases (f c i r > 0.999) 

19 for deriving the clear sky formula and 5788 overcast cases { f c i r < 0.01) for the cloudy sky 

20 formula. These data consist of 43% of all collocated measurements in the training data 

2 1 set. 

22 

23 3.4. Validation 
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1 The modified algorithm has been applied to the CERES Terra Edition 2B and Aqua 

2 Edition 2 A SSF data products. The Terra Edition 2B data, as mentioned previously, spans 

3 from the period from March 2000 to December 2004. The Aqua Edition 2A data covers 

4 the period from July 2002 to March 2005. Fig. 8 illustrates the comparison of the SDLW 

5 of Terra Edition 2B computed with the modified Zhou-Cess algorithm versus the ground 

6 measured SDLW over 29 sites around the globe stratified for different scene types. The 

7 modified algorithm has smaller bias (less than 1 Win 2 ) than the original algorithm for 

8 most of the scene types except for the coastal scenes where the systematic bias is 7.3 

9 Win" . The scatter plots in Fig. 9 of the derived versus measured SDLW shows that the 

10 systematic bias from water cloud and ice cloud are both reduced and the difference 

1 1 between them is also reduced as compared to those from the original algorithm shown in 

12 Fig. 6. 

13 The modified algorithm (named LW model C or LWC here) has also been applied 

14 to Aqua processing along with LWA and LWB, with the following results covering the 

15 33 months from July 2002 to March 2005. The ground measurements are collected from 

16 the same 29 global sites as shown in Fig. 4. All data available from these periods are used 

17 except for some data gaps. The validation was carried for clear sky and cloudy sky 

18 separately. 

19 Table 1 compares the bias and random error of the original algorithm to those of 

20 the modified algorithm as well as LWA and LWB for clear sky fluxes. Using the 

21 modified algorithm reduces the huge negative bias for the ( Antarctic region from -137.0 

22 Win’ to an acceptable -10.5 Win' . The use of the modified algorithm has also resulted in 

23 improvement to either the bias or random error or both for most of the regions. The 
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1 modified algorithm does produce biases for coastal and (Ant)arctic areas of the order of 

2 10 Wm' , whereas globally, there is a less than 1 Win bias and an 18 Wm" random 

3 error. Comparison of LWC with LWA and LWB for clear sky fluxes shows that the 

4 modified LWC has smaller bias for continental region than both LWA and LWB. LWC 

5 does, however, have larger systematic bias in the coastal and desert areas. The global 

6 mean bias is smaller for LWC probably due to cancellation of positive and negative 

7 biases. The random error of 18.5 Wm" is slightly higher than those of LWA (17. 7 Wm" ) 

8 and LWB (16.8 Win" 2 ). 

9 Table 2 compares the bias and random error of the original algorithm to those of 

10 the modified algorithm as well as LWB for cloudy sky fluxes (LWA only computes clear 

11 sky flux). The results in Table 2 show that the major improvement for cloudy sky flux 

12 also occurs at (Ant)arctic region. The bias is reduced from -46.6 Wm" to -0.89 Wm" - . 

13 There is also a modest improvement for the continental areas, though the bias is slightly 

14 larger for coastal and desert regions. Globally, the bias is 1.2 Wm" with a random error 

15 of 23.6 Wm' . Comparing LWC cloud sky flux with LWB, we find that LWC performs 

16 slightly better in continental, (Ant)arctic, desert and island regions, while only slightly 

17 worse in coastal areas. 

18 Similar results for all sky fluxes (clear plus cloudy cases) are shown in Table 3. 

19 The modified algorithm shows improvements over original algorithm for most 

20 geographical regions in systematic error, except for coastal and island regions. 

21 Meanwhile, the random error is reduced for every region. The modified LWC has slightly 

22 smaller systematic bias than LWB for all the regions except for the coastal region, which 

23 still has a systematic bias of 6.6 Wm" . The global mean bias of all sky fluxes is 1.1 Win'" 
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1 and the random error is 23.0 Wm' 2 , slightly better than LWB. Although LWC and LWB 

2 are very different algorithms, we have observed that the random errors of both models are 

3 fairly close to each other for each geographical region. This might be related with spatial 

4 temporal collocation of satellite measurements with ground measurements and the 

5 uncertainties associated with satellite retrieved cloud parameters, i.e., cloud fraction used 

6 in both algorithms. 

7 Despite the overall good performance, there are a few regions or stations that still 

8 present large biases. It is found that the large bias of coastal area mainly comes from the 

9 COVE (Chesapeake light house, 36.90° N, 75.71° W) site which has a positive bias of 

10 18.6 Win 2 . Besides the COVE site, a few other sites, namely, Fort Peck (FPK, 48.31° N, 

1 1 105.10° W), Sioux Falls (SXF, 43.73° N, 96.92° W), Payeme (PAY, 46.82° N, 6.95° E) 

12 and the South Pole (SPL) have systematic bias larger than 10 Wm-2. The causes for the 

13 large biases in these sites maybe related with specific conditions of the sites and should 

14 be examined individually. 

15 Gupta et al. [2004] found that systematic bias of LWB varies with instrument network 

16 and day-night skin temperature bias. To rule out that some of the systematic bias might 

17 be due to instrument setup, we computed the biases grouped with different network. 

18 Table 4 shows no discernible trends in systematic bias between different networks. The 

19 ARM instruments have smaller random error because four out of five instruments are 

20 clustered in the ARM SGP site. The BSRN instruments have the largest random error due 

21 to varied geographical locations of the measurements. The small systematic bias should 

22 not be contributed to instruments themselves. 
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1 The systematic differences between daytime and nighttime biases, however, are found 

2 for all sky fluxes from both Terra and Aqua platforms (Table 5). The revised algorithm 

3 overestimates the daytime SDLW by 3-4 Wm' and underestimates nighttime SDLW by 

4 2-3 Win' . This might be related to the use of skin temperature estimates. The original 

5 algorithm uses 2-meter air temperature which is more closely related to the lower 

6 atmosphere and hence the SDLW. An investigation is currently underway at NASA 

7 Langley Research Center to use a constrained surface temperature in the CERES surface 

8 radiation product since all of the LW models tend to overestimate the surface downward 

9 LW radiation for cases where the surface temperatures are significantly higher than the 

10 air temperatures. The modified skin temperature (or constrained temperature) is expected 

11 to reduce the day-night bias and also the overall bias in this algorithm and others. 

12 

13 4. Conclusions 

14 An improved version of Zhou-Cess algorithm has been formulated which avoids the 

15 large errors in the SDLW at low water vapor conditions by adding an offset to the 

16 logarithmic water vapor term. The new algorithm also utilizes cloud fraction and cloud 

17 liquid and ice water paths available from the CERES SSF product to separately compute 

18 the clear and cloudy portions of the flux. The new algorithm has been validated for the 

19 Terra and Aqua satellites against surface measurements at 29 stations around the globe. 

20 The results show significant improvement over the original version and are comparable 

21 or slightly better than the more sophisticated algorithms currently implemented in the 

22 CERES processing. This revised version of Zhou-Cess algorithm will be incorporated 

23 into the CERES operational processing. 
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1 


2 Table 1. Error statistics of longwave clear sky fluxes computed from LWC (original 

3 LWC, modified LWC), LWA and LWB against surface measurements for Aqua satellite 

4 from July 2002-March 2005. 


Sites 

# of Points 


Original LWC 

LWC 

LWA 

LWB 

Continental 

5012 

Bias Wm" 2 (%) 
a Wm" 2 (%) 

6.71 ( 2.38) 
20.9 ( 7.4) 

-0.16(-0.06) 
16.6 ( 5.9) 

-4.48( -1 .59) 
15.8 ( 5.6) 

-6.93(-2.46) 
15.4 ( 5.5) 

Coastal 

609 

Bias Wm' 2 (%) 
a Wm" 2 (%) 

16.66 ( 5.84) 
17.8 ( 6.3) 

10.03(3.52) 
14.8 ( 5.2) 

4.92 ( 1.72) 
12.9 ( 4.5) 

-0.15(-0.05) 
13.2 ( 4.6) 

Ant(arctic) 

903 

Bias Wm" 2 (%) 
a Wm" 2 (%) 

-137.0(122.0) 
34.5 ( 30.7) 

10.45(9.30) 
11.1 ( 9.9) 

-16.0(14.30) 
11.0 ( 9.9) 

-8.83( -7.88) 
11.0 ( 9.9) 

Desert 

1640 

Bias Wm" 2 (%) 
a Wm" 2 (%) 

-1.68 ( -0.53) 
24.4 ( 7.8) 

-7.12(-2.27) 
23.9 ( 7.6) 

-0.41 (-0.13) 
22.8 ( 7.3) 

-5.15 (-1.64) 
21.1 ( 6.7) 

Island 

138 

Bias Wm" 2 (%) 
a Wm' 2 (%) 

9.31 ( 2.47) 
10.3(2.7) 

2.92( 0.77) 
12.6 ( 3.3) 

-0.75 (-0.20) 
12.0 ( 3.2) 

-0.77 (-0.20) 
13.8 ( 3.6) 

Global 

8302 

Bias Wm' 2 (%) 
a Wm" 2 (%) 

-9.81 (-3.61) 
41.4 (15.3) 

0.42 ( 0.15) 
18.5 ( 6.8) 

-4.18 (-1.54) 
17.7 ( 6.5) 

-6.19 (-2.28) 
16.8 ( 6.2) 


5 

6 
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1 


2 Table 2. Error statistics of long wave cloudy sky flux computed using LWC (original, 

3 modified) and LWB against surface measurements for Aqua satellite from July 2002- 

4 March 2005. 


Sites 

# of Points 


Original LWC 

LWC 

LWB 

Continental 

21117 

Bias Wm 7 (%) 
ct Wm" 2 (%) 

-7.42 ( -2.33) 
28.0 (8.8) 

0.61 ( 0.19) 
22.6 ( 7.1) 

-3.38 (-1.06) 
22.4 ( 7.0) 

Coastal 

3302 

Bias Wm' 2 (%) 
a Wm" 2 (%) 

1.80 ( 0.52) 
21.2 ( 6.1) 

6.00 ( 1.72) 
18.5 ( 5.3) 

2.51 ( 0.72) 
19.0 ( 5.4) 

Ant(arctic) 

17663 

Bias Wm" 2 (%) 
a Wm" 2 (%) 

-46.63 (-20.17) 
50.6 (21.9) 

-0.89 (-0.39) 
23.6 ( 10.2) 

-6.17 (-2.67) 
25.0 ( 10.8) 

Desert 

4170 

Bias Wm" 2 (%) 
a Wm" 2 (%) 

8.19 ( 2.42) 
28.1 (8.3) 

9.67 ( 2.85) 
27.7 ( 8.2) 

11.88 ( 3.51) 
28.8 ( 8.5) 

Island 

6729 

Bias Wm" 2 (%) 
ct Wm' 2 (%) 

0.58 (0.14) 
13.5 (3.3) 

0.95 ( 0.23) 
12.5 ( 3.0) 

5.71 ( 1.39) 
14.8 ( 3.6) 

Global 

52981 

Bias Wm' 2 (%) 
ct Wm" 2 (%) 

-17.67 (-5.80) 
40.1 (13.2) 

1.20 ( 0.39) 
23.6 ( 7.8) 

-1.59 (-0.52) 
24.4 ( 8.0) 


5 

6 
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1 


2 Table 3. Error statistics of long wave all sky flux computed using LWC (original, 

3 modified) and LWB against surface measurements for Aqua satellite from July 2002- 

4 March 2005. 


Sites 

# of Points 


Original LWC 

LWC 

LWB 

Continental 

26129 

Bias Wm 7 (%) 
a Wm" 2 (%) 

-4.71 (-1.51) 
27.2 ( 8.7) 

0.46 ( 0.15) 
21.6 ( 6.9) 

-4.06 (-1.31) 
21.3 ( 6.8) 

Coastal 

3911 

Bias Wm' 2 (%) 
a Wm" 2 (%) 

4.12 ( 1.21) 
20.9 ( 6.1) 

6.63 ( 1.95) 
18.1 ( 5.3) 

2.10 ( 0.62) 
18.4 ( 5.4) 

Ant(arctic) 

18566 

Bias Wm" 2 (%) 
o Wm" 2 (%) 

-51.03 (-22.64) 
51.6 (22.9) 

-0.34 (-0.15) 
23.3 ( 10.3) 

-6.30 ( -2.79) 
24.8 (11.0) 

Desert 

5810 

Bias Wm" 2 (%) 
a Wm" 2 (%) 

5.40 ( 1.63) 
27.5 ( 8.3) 

4.93 ( 1.48) 
27.4 ( 8.3) 

7.08 ( 2.13) 
27.3 ( 8.2) 

Island 

6867 

Bias Wm" 2 (%) 
a Wm' 2 (%) 

0.76 (0.19) 
13.5 (3.3) 

0.99 ( 0.24) 
12.5 ( 3.1) 

5.58 ( 1.36) 
14.8 ( 3.6) 

Global 

61283 

Bias Wm' 2 (%) 
a Wm" 2 (%) 

-16.61 (-5.54) 
41.1 ( 13.7) 

1.09 ( 0.36) 
23.0 ( 7.7) 

-2.22 ( -0.74) 
23.5 ( 7.8) 


5 
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1 

2 Table 4. Bias of long wave all sky fluxes computed using LWC compared with surface 

3 measurements stratified with observational network. 


Platfonn 


ARM 

SURF 

BSRN 

CMDL 

Terra 

Bias Wm 2 

1.1 

2.9 

-3.1 

0.5 


a Wm" 2 

17.6 

22.2 

26.2 

21.2 

Aqua 

Bias Win' 2 

2.7 

0.9 

-0.2 

-1.1 


a Wm" 2 

18.8 

22.6 

27.6 

21.1 


4 

5 
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1 


2 Table 5. Bias of long wave all sky fluxes computed using LWC compared with 

3 surface measurements stratified with observational time. 




Daytime 

Nighttime 

Terra 

Bias Wm 2 

4.9 

-3.4 


a Wm" 2 

15.7 

24.3 

Aqua 

Bias Wm 2 

3.5 

-2.3 


a Wm' 2 

20.9 

24.2 


4 

5 
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1 




3 Fig. 2 Relationships between MODTRAN computed clear sky SDLW and SULW (upper 

4 panel) and flux ratio versus column water vapor (lower panel). The black and dot lines 

5 are linear regression fits (upper panel) and log-square fits (lower panel) using data points 

6 with PWV greater than 0.5 cm (solid line) and PWV less than 0.5 cm (dot lines), 

7 respectively. 
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1 

2 Fig. 3 Normalized probability distribution functions of SDLW, SULW and PWV from 

3 global network (solid line), ARM SGP site (dash-dot line) and ARM NSA site (dot line). 
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Fig. 4 Ground stations used in this study. 
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3 Fig. 5 Scatter plots of SDLW computed using original Zhou-Cess algorithm for Terra 

4 satellite versus ground measurements for different scene types. 
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2 Fig. 6 Scatter plots of SDLW computed using original Zhou-Cess algorithm versus 

3 ground measurements stratified for clear sky (Area_clr > 99.9%), water cloud (LWP>5 

4 g/nr) and ice cloud (LWP < 5 g/nr, IWP>1 g/nr) conditions for Terra satellite. 
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3 Fig. 7 Logarithmic normalized probability distribution functions of cloud liquid water 

4 path from ARM SGP site as observed from ARM microwave radiometer (dotted line) and 

5 CERES instruments (solid line). 
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2 

3 Fig. 8 Scatter plots of SDLW computed using modified Zhou-Cess algorithm versus 

4 ground measurements for different scene types for Terra satellite from March 2000 to 

5 December 2004. 
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1 

2 Fig. 9 Scatter plots of SDLW computed using modified Zhou-Cess algorithm versus 

3 ground measurements stratified for clear sky (Area_clr > 99.9%), water cloud (LWP>5 

4 g/m 2 ) and ice cloud (LWP < 5 g/m 2 , IWP>1 g/m 2 ) conditions for Terra satellite from 

5 March 2000 to December 2004. 
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